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Thermal response functions of strongly correlated electron systems are of appreciable interest to 
the larger scientific community both theoretically and technologically. Here we focus on the infinitely 
correlated t-J model on a geometrically frustrated two-dimensional triangular lattice. Using exact 
diagonalization on a finite sized system we calculate the dynamical thermal response functions 
in order to determine the thermopower, Lorenz number, and dimensionless figure of merit. The 
dynamical thermal response functions is compared to the infinite frequency limit and shown to be 
very weak functions of freq uency , hence, establishing the validity of the high frequency formalism 
recently proposed by ShastryH^El for the thermopower, Lorenz number, and the dimensionless figure 
of merit. Further, the thermopower is demonstrated to have a low to mid temperature enhancement 
when the sign of the hopping parameter t is switched from positive to negative for the geometrically 
frustrated lattic^ considered. 

PACS numbers: 72.15.Jf, 65.90.+i, 71.27.+a 



I. INTRODUCTION 

There is current interest in the physics, as well as in- 
dustrial and engineering communities regarding thermo- 
electrics of strongly correlated electron systems. This in- 
terest has been recently revived by the demonstration of 
the unexpectedly high thermopower s een i n the very in- 
teresting material sodium cobalt oxid^^'^'^NCO). Theo- 
retically, thermoelectrics have been a long standing prob- 
lem in physics especially when concerned with strongly 
correlated systems which are not amenable to perturba- 
tive treatments. 

Essentially there are two standard theoretical ap- 
proaches concerning the problem of thermoelectrics, cf. 
Ref. I8I9I10I1I1 The first makes use of Boltzmann theory 
often complimented with standard Fermi liquid theory. 
This methodology is reliable for weakly coupled prob- 
lems where long lived quasiparticles remain well defined 
and where perturbation theory remains valid. The sec- 
ond approach is to use the full rigor of the Kubo formal- 
ism which is valid for all situations but whose dynamical 
character make it unyielding and difficult to make any 
real progress, especially for strongly correlated systems. 

Recently, ShastrjEEIll has proposed a third method 
which handles the strong electron interactions with the 
respect which they deserve while avoiding the complexity 
of the full dynamics of the Kubo formalism. Very briefly, 
this methodology considers the ultimate quantities one 
is usually interested in when calculating conductivities. 
Often there is more interest in certain combinations of 
conductivities, which form more experimentally accessi- 
ble quantities (such as the thermopower or Seebeck coeffi- 
cient, Lorenz number, dimensionless figure of merit, the 
Hall coefficient, etc.), than in the conductivities (elec- 
trical, thermal, etc.) themselves. The basic proposal 
is that certain combinations of conductivities have weak 
dynamical character (weak frequency dependence), and 
thus lend themselves to a high frequency expansion. The 
upshot of this expansion is that it yields formulas that 



are much simpler, although non-trivial, than the Kubo 
formulas and yet the interactions are fully respected com- 
pared to the usual approximations which risk missing im- 
portant effects. 

Nearly 15 years previously, this basic high frequency 
expansion methodology was originally employed by Shas- 
try, Shraiman, and SinghP^ to calculate the Hall coef- 
ficient (at high temperatures) for a strongly correlated 
electron model {t-J mo del) with success. In the last 
year, the present authors^'^^-^' have applied this high fre- 
quency expansion to calculate the Hall coefficient and 
thermopower for the very interesting NCO system ex- 
plaining in both quantitative and qualitative detail the 
physics of the so-called Curie- Weiss metal for Naa;Co02 
at electron doping x ~ 0.7. This Curie- Weiss metal dis- 
plays behavior that is an interesting hybrid between those 
of insulating and metallic systems. The high frequency 
formalism allowed the investigation of this complicated 
system successfully by incorporating the important ef- 
fects of interactions. 

Interestingly, Shastry, via the high frequency expan- 
sion, was able to predict a low to mid temperature ther- 
mopower enhancement due to a change in sign of the 
hopping parameter t of the t-J model for a geometrically 
frustrated^ two-dimensional triangular lattic4^. This 
corresponds to a fiduciary hole doped C0O2 layer of NCO 
that has yet to be realized experimentally. The lattice 
topology plays a crucial role in this enhancement as it 
owes its existence primarily to electron transport and is 
not thermodynamic or entropic in origin. This predic- 
tion was put on firmer footing by the present authors 
in Ref. [THl concerning the Curie- Weiss metallic phase of 
NCO which itself has an underlying two-dimensional tri- 
angular lattice. We emphasize that we work with a hole 
doped system (0 < n < 1, n electron density) and in 
order to compare with experiments on NCO, w e perform 
a suitable particle-hole transformatioEp'^IISIli! 

In this work, we establish the validity and accuracy 
of the high frequency formalism for the thermopower. 
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Lorenz number, and figure of merit for the strongly cor- 
related electron t-J model on a two-dimensional triangu- 
lar lattice. This is accomplished by comparing the high 
frequency expressions with those obtained via the full 
Kubo formalism. This comparison for a strongly corre- 
lated system, such as our model, is only possible through 
numerical exact diagonalization of a relatively small sys- 
tem (£ = 12 site lattice). However, we feel that our re- 
sults provide a much desired and important benchmark 
for the high frequency formalism establishing its effec- 
tiveness and usefulness. 

Furthermore, the t-J model is generally representa- 
tive of strongly correlated electron models so our results 
should be applicable to other strong correlation models 
such as the Hubbard model for large U. While the geo- 
metrically frustrated triangular lattice provides an inter- 
esting enhancement of the thermopower for intermediate 
temperatures, the general validity of our results should 
obtain for other lattice topologies (frustrated and non- 
frustrated) . 

The plan of this paper is as follows: In section [TT| we 
describe the details of our model and the exact diago- 
nalization used, section |HI| generally quotes the Kubo 
formulas for the considered conductivities and the hi gh 
frequency formulas from ShastrjW2l. In sections rv|vT 



report results for the thermopower (for both positive and 
negative hopping t), the Lorenz number, and the figure 
of merit (both for positive hopping t > 0), respectively. 
Section |VII| concludes while some formulas are given for 
completeness in the appendix. 



II. t-J MODEL AND DIAGONALIZATION 

As mentioned above we study the t-J model Hamil- 
tonian which describes a strongly correlated hole doped 
Mott insulator. The Hamiltonian is 
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Jirf)Sp- Sp+ff , (1) 



where ct^(cfo-) — PGCpai^fcr)PG are Gutzwiller projected 
fermion creation(destruction) operators where the pro- 
jection operator Pq projects out all doubly occupied lat- 
tice sites. The lattice vector rf connects nearest neighbors 
which are coupled (with strength J{f])) via their spin de- 
gree of freedom (Sp is the three-component spin opera- 
tor). For simplicity we take the hopping t{ff) = t and 
spin J(r/) = J coupling parameters to be constants and 
the lattice constant has been set to unity. 

In this work we also aim to apply our calculations 
to the experimental system of NCO which is electron- 
doped andhas been previously modeled using the t-J 
j-^Q(jg{i3|i5|i6li71 symmetry of the Hubbard 

model with regard to half filling to map our system to 
NCO, i.e., we apply the replacement rules t — > —t, dop- 
ing X — |1 — n|, and —q^ where n is electron density 
per site and q^ = —\e\ is the electron charge. At this 




FIG. 1: C = 12 site torus geometry used throughout this 
work. We diagonahze the Hamiltonian at all densities of this 
lattice for the t-J model, however, as explained in the text we 
label the electron densities in terms of hole doping x. 



point we will abandon referencing particular systems by 
the electron density and instead reference them by the 
doping X = n ~ I. 

Thermodynamics will be considered within the canon- 
ical ensemble and considering that the goal of this paper 
is to calculate full thermodynamic Kubo formulas it is 
a prerequisite that we obtain the full eigenspectrum of 
the system in order to calculate all relevant current ma- 
trix elements. Therefore we make progress through the 
exact numerical diagonalization of a finite system. At 
this level of study various approximations such as the 
finite temperature Lanczos methocP^, dynamical mean 
field theory^, etc., could perhaps muddy the issue of 
comparing the high frequency expansion of various ther- 
moelectric properties to their full Kubo formulations and, 
hence, will not be considered. Further, we are interested 
in the behavior of the system over a whole range of densi- 
ties and do not wish to confuse the results by considering 
a smattering of small system sizes along with their partic- 
ular inherent and unavoidable finite size effects. There- 
fore, we find that the largest two-dimensional lattice that 
we can fully diagonalize in a satisfactory way is a £ = 12 
site toroidal lattice, c.f. Fig[T] Another reason for using 
this lattice is that it w as use d extensively in the present 
authors' previous work^^^'^^ on the two-dimensional tri- 
angular lattice t-J model applied to NCO. 

To reduce the computational demand of the exact diag- 
onalization to a more manageable size we employ a num- 
ber of symmetries. Firstly, we consider only the largest 
^2 sector of the full Hilbert space, i.e., the smallest \Sz\ 
subspace which is \Sz\ — 0{Ti/2) for even(odd) number 
of electrons. This sector of the full Hilbert space domi- 
nates the physics so this "approximation" is as good as 
exaciP^. The most useful symmetry employed is trans- 
lational invariance which essentially reduces the Hilbert 
space dimension by a factor of C. Our largest Hilbert 
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space dimension occurs for x = 1/3 corresponding to 
34650 states in the |S'^| = subspace and after applying 
translational symmetry we need only diagonalize matri- 
ces of dimension ~ 2900. While this matrix dimension is 
not particularly huge in the realm of matrix diagonaliza- 
tion we must consider a double sum over these ^ 2900 
states to calculate the Kubo formulas. This double sum 
is quite time consuming and limits our abilities to con- 
sider larger lattices and, in fact, it limits our abilities to 
consider all doping values x on the chosen £ = 12 site 
lattice used here. The dopings where we can calculate 
the full Kubo formula is limited to x > 0.5 and x < 0.2. 



III. DYNAMICAL THERMAL RESPONSE 
FUNCTIONS 

In this section we quote the formulas for the Kubo 
linear response for thermoelectrics following very closely 
the work of Shastrjffl^El particular we are interested in 
the electrical a{uj,T), thermoelectrical 7(0;, T), and the 



thermal k{uj,T) conductivities, respectively. In terms of 
these conductivities the ultimate goal of this work are 
commonly measured physical quantities of interest such 
as the thermopower S, the Lorenz number L, and the 
dimensionless figure of merit ZT commonly defined^^i^^ 
as 



and 



(2) 



L(c.,T) = ;^(^-{5(c.,T)}^ (3) 
{S{u;,T)r 



Z{lo,T)T 



L{uj,T) 



(4) 



These conductivities all have familiar Kubo formulas 
which, in the Lehmann representation, are written as 



cr(wc,T) 



huJrC 



n \ ^ e 



huj, 
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flUJr 



{n\J^\m){m\J^\n) 
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and 



K,{u)c,T) 



hujcTC 



Z £n - £,n + Tt-I^c 



|HjQ|m)|^ 



(7) 



In the above, \k) is a normalized eigenstate of the Hamil- 
tonian with energy e^, Z — exp(— /3ej.) is the canon- 
ical partition function, and (3 = l/ksT is the inverse 
temperature. A thermal average is indicated by (•••). 
The dynamical temperature variation is turned on adia- 
batically from the infinite past, i.e., tOc = uj + iO"*". 

In Eqs. [Hj [rjthe charge current is formally given by 

Jx^~ lini -l^[K{kx),qen{-kx)] (8) 

while the heat current J§ is 

JQ = -lhn^l^jK{kx),k{-kx)] . (9) 

Here K = H—^ih is the grand canonical Hamiltonian and 
^(^) = X)^^P(*^ ■ r)^(r) is the Fourier decomposition 
for mode A: of a local operator A{r) which operates at 



position r. It should be noted that the heat and charge 
currents are related via the energy current as — — 
{fi/qe)Jx with 

if = - ^hm^ ^jT{kx), ln~kx) + U{-kx)] , (10) 

with T and U being equal to the kinetic and potential 
energy operator^^U^ respectively. 

From Ref. lll2l51 the definitions for stress tensor f^x, the 
thermoelectric operator ^xx, and the thermal operator 
Qxx are 

Txx = - lim -^[Jx{kx),qen{-kx)] (H) 

^xx = - ^[Jx{kx),K{-kx)] (12) 
k^^a akx 

Qxx = - lim -^[jQ{kx),k{-kx)] . (13) 
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The explicit form of these operators for the t-J model 
are given in the appendix. 

Throughout this work we take the chemical potential 
/i(T) in the above formulas to be that obtained from the 
canonical ensemble, i.e., /i(T) = dF/dN, where F = 
— (l//3)ln(Z) is the Helmholtz free energy. For a finite 
sized system we approximate this partial derivative as 
fi{T) = {Fn+1 - Fn-i)/2 for an N electron systeniSD. 

When calculating the full frequency dependent conduc- 
tivities for finite sized clusters one must take into account 
the discreteness of the energy spectrum caused by the fi- 
nite nature of the cluster. This is done by introducing 
a broadening factor 77 where the frequency then becomes 
LDc ^ Lo -\- irj. The broadening factor is taken to be the 
mean energy spacing between states with non-zero cur- 
rent matrix elements. Table [l] provides values of rj for 
the systems considered here and is generally weakly 
dependent on x and of order 77 ~ 3|t|. 



TABLE I: Broadening factor rj for the t-J model on a two- 
dimensional triangular lattice with C = 12 sites. Both pos- 
itive and negative values of the hopping t are given. The 
value of J has very little effect on 77 compared to the weak x 
dependence. 



0.83 
0.75 
0.67 
0.58 
0.17 
0.083 



■n/\t\ (J = Q.2\t\,t > 0) 



3.229589 
4.198734 
4.659046 
4.922330 
3.786627 
2.803789 



■n/\t\{j ^ Q,t < 0) 



2.678842 
4.221241 
4.697231 
4.934004 
3.772738 
2.769104 



Following the work of Ref. I1|2|3I one can consider 
the high frequency limit of the thermopower, Lorenz 
number, and figure of merit in the hope and expec- 
tation that for strongly correlated systems modeled by 
the t-J model these combinations of conductivities will 
have weak frequency dependence and yet still capture 
the essential strongly correlated physics, i.e., the Mott- 
Hubbard physics. For notational convenience we indicate 
the high frequency expansions by a superscript star, that 
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S*{T) = lim S{uj,T) = 



(14) 



L*{T) = lim L{uj,T) = 



-{S*{T)r, (15) 



and 



Z*{T)T^ lim Z{uj,T)T = 



{s*{T)r 

L*{T) 



(16) 



While the high frequency starred quantities are not 
trivial to calculate they are considerably simpler than 
the full dynamical Kubo formulas (Eqs. 2]|4| as they are 



equilibrium expectation values and not dynamical in na- 
ture. It is also reasonable to expect that in the future, 
approximate methods could be used to calculate their 
full temperature dependence that are not as limited as 
the exact diagonalization brute force methods used here. 
However, one must be very careful when aiming to estab- 
lish new approximation for strongly correlated systems 
and, hence, the need for the brute force calculations in 
this work. 



IV. THERMOPOWER 

The thermopower can be factored instructiveljff^ as 
7(w,T) 7(w,r) ii.{T) 



S{uj,T) 



(17) 



where -7(10, T) = 7(w,r) - {^{T) / q^T)(j{uj,T) defines 
7(w,r). This factorization displays clearly the two con- 
tributions composing the thermopower; one arising from 
electron transport and the other from thermodynamics 
(entropy). This is discussed in detail below. 

As was mentioned previously, calculating this type of 
formula for strongly correlated systems is very difficult. 
Often there is a desire to give the transport term little 
importance and, therefore, drop it. This leaves merely 
the second term in Eq. [l7] involving the chemical poten- 
tial alone. This is the so-called Mott-Heikes (MH) term 
for the thermopower which is valid at high temperatures 
and is described in more detail below. However, as will 
be shown below, for low to intermediate temperatures 
this approximation is not adequate. 

On physical, as well as theoreticaP^, grounds, one ex- 
pects the thermopower to vanish at T = and for a non- 
interacting system it is simple to show this. For our pur- 
poses it is instructive to describe this vanishing through a 
delicate balancing act where the transport term exactly 
equals the MH term as the temperature tends towards 
zero, i.e.. 



lim .7(^„Mm\ 







(18) 



Achieving this balance in a finite sized system is not 
possible explicitly, even for non-interacting electrons. 
However, it does suggests a formulation of the ther- 
mopower into the two contributing terms mentioned 
above; a frequency dependent transport term Str{oJ,T) 
and a frequency independent Mott-Heikes (MH) term 
Smh{T) with both defined through 



5(a;,T) = 



1 \ T^{u:,T) T7(c^,0) 
T\a{uj,T) a(c^,0) 

f /i(r)-M(0) ] 

I 9eT J 

Str{0J,T) + SMH{T) . 



(19) 



Therefore, even for a finite sized system one can obtain 
a transport term and MH term that independently equal 
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zero at T = ensuring that S{u!,T) vanishes in the zero 
temperature Hmit as expected. This is not the complete 
picture however. The MH term contains the chemical 
potential which is expected to behave quadratically in T 
as T — > for thermodynamically large systems. Finite 
systems have two particular differences. One is that the 
spectrum is discrete giving rise to a ground state energy 
gap in situations without degeneracies. In those instances 
there will be a low temperature exponential behavior of 
/i(T) which is not really a problem for our purposes be- 
cause fi(T)/T will still vanish in the zero temperature 
limit. The existence of ground state degeneracies, on the 
other hand, are a bigger concern. Their existence cause 
the chemical potential to behave linearly in T at low tem- 
peratures which, in turn, produces a MH term that does 
not vanish. We argue that this is an unwanted unphysical 
result and our solution is to merely discount this ground 
state degeneracy (when it exists) when calculating fi(T) 
ensuring that the MH term vanishes as T ^ 0. 

The high frequency expansion of S{lu,T) is similarly 
written as 



S*{T) 



Sl{T) + SMHiT). 



+ SMHiT) 



(20) 



Again we have defined similarly to 7(0;, T) through 

^xx *^xx if^i'-^^ / Qe)^xx ■ 

The transport term of the thermopower eventually 
vanishes as T becomes large so we know that the Mott- 
Heikes term eventually dominates the thermopower and 
becomes useful for a number of reasons. One reason it is 
that it is not a dynamical quantity and, hence, is easier to 
compute. Secondly, it is often hoped that the MH term 
dominates the thermopower and one only needs to con- 
sider it. This is due to the fact that at high temperatures 
it approaches a constant since /^(T) is eventually linear 
in T. Previous work by BenP^, and Chaikin and BenP^ 
worked out the infinite temperature limit of Smh{T) for 
a number of systems. There is an elegance and simplicity 
to these formulas since the infinite temperature limit of 
Smh{T) is determined merely from counting arguments 
related to the Hilbert space dimension of the problem at 
hand. A central question regarding the MH term is how 
low of a temperature does the MH limit remain a valid 
approximation to the full thermopower. We provide an 
answer to that question for the t-J model in this work 
which is discussed later. 

The two MH limits we consider he re are f or the uncor- 
related band^^ and for the t-J modeyJ^SllZI 



lim Smh{T) 



In 

Inf^ 
a In [ 



n 

2-ri 



uncorrelated with < n < 2 
t-J with < n < 1 
t-J with 1 < n < 2 



(21) 



remembering that 



is the electric charge. 



Although a somewhat blunt formulation, the Mott- 
Heikes limits already contain a plethora of information. 
For example, for the uncorrelated model, even with a fi- 
nite interaction parameter U, the thermopower diverges 
at X ^ 1 (n — > 0), is positive for all hole dopings a; > 
(n < 1) and is exactly zero for the half filled case a; = 
(n = 1). For electron doping the thermopower would be 
purely negative diverging negatively as x ^ 1 (n — > 2). 
In the whole range of densities (0 < n < 2) the MH limit 
would predict a single sign change. 

For the t-J model the essentially infinite strength in- 
teractions of the electrons causes two additional sign 
changes compared to the uncorrelated or finite U Hub- 
bard model. The thermopower still diverges as a; ^ 1, 
is positive for x > 1/3, negative for x < 1/3, and di- 
verges negatively at half filling. For electron doping we 
use particle hole symmetry to get precisely the opposite 
behavior: a positive divergence at half filling, positive 
for X < 1/3, a sign change to negative thermopower for 
x > 1/3 and a negative divergence as a; ^ 1. Hence, two 



additional zero crossings emerge due to the interactions. 

In the following subsections we report results for 
S*{T), the difference S'(w,T) - S*{T) and the Mott- 



Heikes term S mh(T ) for both positive (Sec. IV 1 1 and 
negative (Sec. IV 2 1 signs of the hopping t. In all fig- 



ures the thermopower is given in experimental units of 
^iV/K where ks/lqel — 86fiV/K, and we have multiplied 
by (—1) to facilitate comparison with the electron doped 
NCO system, see Eq. [21] 



1. Positive hopping t > 

Fig. [2] shows S*{T) as a function of both doping and 
temperature for the case of positive hopping t > 0. 
We have computed S*{T) for two different values of J, 
namely, J = (Fig. [2ji) and J = 0.2\t\ (Fig. |2|d). Pro- 
jected onto the T — 10\t\ plane are the two MH limits; 
the uncorrelated model (blue dashed line) and the t-J 
model (solid red line) . The high temperature behavior of 
S*{T) matches the t-J model MH limit expectation quite 
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t>0, J=Oltl 



s'(|lV/K] 



(a) 




(b) 

FIG. 2: (color online) S*{T) as a function of doping x and 
temperature T for positive hopping t > corresponding to 
NCO after particle-hole transformation. Panel (a) and (b) 
are for J — and 0.2|t|, respectively. Projected onto the 
T = 10|i| plane is the Mott-Heikes limits for the uncorrelated 
(blue dotted) and the t-J models (red), respectively. S*{T) 
approaches the MH limit for the t-J model relatively quickly, 
i.e., by approximately T ~ 6|t|. The horizontal black dotted 
lines indicate the position of zero thermopower. 



satisfactorily producing a sign change near x — 1/3. The 
slight difference between S*{T = I0\t\) is a combination 
of T = 10|i| being large but finite and finite size effects 
coming from the finite dimension of the Hilbert space 
yielding a discrepency even at T = oo. 

For large dopings x > 0.58 the thermopower mono- 
tonically grows to a somewhat large value of (100 — 
200)fiV/K, growing faster with temperature the higher 
the doping. The value of J has little to no effect in this 
range of doping. The thermopower is pinned at zero for 
T = and needs to eventually increase to its MH limit 
which in this range of doping is positive and large. Since 
the doping is large there is little interaction between elec- 
trons and evidently they effectively avoid one another. 
Hence the transport term S'jj,(T) has very little impact 



t>0, J=Oltl 



(a) 




(b) 

FIG. 3: (color online) The Mott-Heikes term of the ther- 
mopower SMniT) versus doping x and temperature T for 
positive hopping t > corresponding to NCO after particle- 
hole transformation. Panel (a) and (b) are for J = and 
0.2jfj, respectively. The black line is Smh{T) while the red 
line is the full S*{T). The MH limits are projected onto the 
T = 10\t\ plane. 



on the full thermopower. This physics is borne out by 
comparing S*{T) to only the Mott-Heikes term Smh{T) 
which is shown in Fig. |3]i-b where there is very little dif- 
ference between S*{T) and Smh{T). 

For dopings x < 0.5 strong electron correlation ef- 
fects are obtained. This is due to the transport term 
which acts to reduce the thermopower at low to inter- 
mediate temperatures. S*{T) no longer monotonically 
approaches its MH limit. Again this is shown more dis- 
tinctly when one compares S*{T) to Smh{T) in Fig.|3|L-b 
where the MH term overestimates the thermopower in- 
dicative of a very active and important transport term 
5j^(T) which serves to reduce the thermopower. Eventu- 
ally, the transport term vanishes as T becomes large and 
the MH term again dominates. 

As the doping approaches half filling [x 0) the ther- 
mopower begins to be purely negative and nearly mono- 
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x=0.83, t>0, J=0.2ltl 



x=0.75. t>0, J=0.2ltl 



(S((0)-S )[HV/K] 



(S((0)-S )[HV/K] 





(a) 



(b) 



x=0.67, t>0, J=0.2ltl 



x=0.58, t>0, J=0.2ltl 



(S((0)-S )[HV/K] 



(S((0)-S )[HV/K] 





(c) 



(d) 



x=0.17, t>0, J=0.2ltl 



x=0.083, t>0, J=0.2ltl 



(S(m)-S )[|1V/K] 




(S(CO)-S )[HV/KJ 

10 r 



0.3 




(e) 



(f) 



FIG. 4: (color online) S{lj,T) — S*{T) as a function of frequency lu and temperature T for J = 0.2|t|, positive sign of the 
hopping t > (corresponding to NCO after particle-hole transformation) and for dopings (a) x — 0.83, (b) x — 0.75, (c) 
X = 0.67, (d) X = 0.58, (e) x — 0.17, and (f) x — 0.083. The frequency dependence is evidently quite weak for dopings x > 0.58 
((a)-(d)). For x =0.17, and 0.083 ((e) and (f)) there is a much stronger frequency dependence that occurs at extremely small 
u! and T. This is most likely due to finite size effects of our £ = 12 site lattice and not an intrinsic property of the t-J model at 
these dopings. For parameters ui > 3\t\ and T > 2\t\ the frequency dependence is approximately flat. The doping of a; = 0.92 
corresponds in our case to only one electron and, hence, has no frequency dependence. 



tonically approaches its now negative MH limit. Some of 
the violent behavior at the lowest temperatures reported 
is no doubt due to peculiarities of the finite sized lattice. 
That aside, the transport term has an increasingly im- 
portant role to play at low to intermediate temperatures 



as the doping is reduced. Again, this is quite obvious in 
Fig.[3^-b where the MH and transport terms are quite di- 
vergent. Hence, it is definitely not a good approximation 
to use SAiniT) as a representative of the thermopower 
in low doping regions of strongly correlated systems for 
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t<0, J=Oltl 




(a) 

t<0, J=0.4ltl 



S [nV/Kl 




(b) 

FIG. 5: (color online) S*{T) as a function of doping x and 
temperature T for negative hopping t < corresponding to a 
fiduciary hole doped C0O2 compound. Panel (a) and (b) are 
for J = and 0.4|t| , respectively. Projected onto the T = 10|t| 
plane is the MH limits. S* (T) approaches the MH limit for the 
t-J model relatively quickly, i.e., by approximately T ~ 6|i|. 
The horizontal black dotted lines indicate the position of zero 
thermopower. 



low to intermediate temperatures (T < 5|t|). 

The value of J has almost no effect on the the ther- 
mopower until a doping of a; = 0.25 is reached. Not 
shown are similar results for J = 0.4|t|. 

It is interesting to note that S*{T) can be well ap- 
proximated by the infinite temperature MH limit for all 
dopings and all J when the temperature is at or above 
approximately 5\t\ < T < 6|i|, as the thermopower has 
an overwhelming MH constribution by that temperature. 

Lastly we justify the use of S*{T) instead of the 
full dynamical thermopower S{uj,T). Fig. |4^-f shows 
S{lu,T) — S*{T) as a function of temperature and fre- 
quency uj for dopings x =0.83, 0.75, 0.67, 0.58, 0.17, and 

0. 08. The dopings x =0.5, 0.42, 0.33, and 0.25 cannot be 
calculated at this time due to computational constraints, 

1. e., the double sum over the current matrix elements in 



t<0, J=Oltl 



(a) 




(b) 

FIG. 6: (color online) The Mott-Heikes term of the ther- 
mopower SMniT) versus doping x and temperature T for 
negative hopping t < corresponding to a fiduciary hole 
doped C0O2 compound. Panel (a) and (b) are for J = 
and 0.4|t|, respectively. The black line is Smh{T) while the 
red line is the full S*{T). The MH limits are projected onto 
the T = 10\t\ plane. 



Eqs. [5]and[6]is quite prohibitive. For x < 0.58 the differ- 
ence between the dynamical thermopower and the infinite 
frequency expansion is less than 2.5fj,V/K and hence has 
very little absolute effect (approximately less than ~ 2% 
difference). Further it should be noted that for values of 
Lu > 3|t| and temperatures T > 2\t\ the frequency depen- 
dence of the thermopower is nearly nonexistent. 

For smaller dopings, i.e., x =0.17, and 0.083, a more 
severe difference between S{u!,T) and S*(T) is found. 
However, this larger difference takes place at extremely 
small frequencies and temperatures especially consider- 
ing this calculation is done for a finite sized system. Re- 
call from Fig. ^ that at temperatures below approxi- 
mately T ~ 0.3|i| the thermopower displays a drastic 
behavior as it approaches T = 0. This extreme behavior 
is almost certainly a consequence of the finite sized lat- 
tice on which we work and not an intrinsic property of 
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the t-J model. In fact, in Fig. |4]3 andj^f we have cut the 
temperature off below T = 0.3|t| as the thermopower is 
badly divergent. Therefore, it should be concluded that 
the frequency dependence of S{uj,T) is most likely weak 
even for very small dopings. 



2. Negative hopping t < 

We now consider the thermopower for a fiduciary sys- 
tem where we have switched the sign of the hopping t, 
i.e., t < 0. Since we study < n < 1, by a particle-hole 
transformation, this corresponds to i > for 1 < n < 2, 
which differs from Na2;Co02 in the sign of the hopping 
(t < for NCO). Such a system does not exist in the lab- 
oratories, and our hope is that our result will stimulate 
the search for a hole doped C0O2 system. 

In view of the topology of the geometrically frustrated 
triangular lattice there is the possibility of important ef- 
fects coming from primarily the transport term. In RefjTl 
Shastry obtained a high temperature expansion of S* (T) 
for the triangular lattice which we quote here for com- 
pleteness for the situation corresponding to < n < 1 
and i > 0, 



S*iT) 



- 1- 



1 



2T 



(22) 



with hole doping x = 1 — n. Comparision with an elec- 
tron doped system such as NCO is facilitated through the 
particle-hole transformation given previously, c.f. Sec.[TTl 
The first term in Eq. [22] is the Mott-Heikes term. The 
second term is due to the transport and serves to reduce 
the thermopower from its large MH upper limit as the 
temperature is decreased. Importantly, it depends on a 
single power of the hopping t. Hence, if one could switch 
the sign of the hopping then S* (T) would evidently grow 
to a maximum as the temperature was decreased before 
approaching its zero value at T = 0. A term similar to 
the second term in Eq. |22] with an odd power of t is only 
present in cases where the topology of the underlying 
lattice is geometrically frustrated. 

Fig. [5)i-b shows S*{T) for the case of negative hop- 
ping as a function of x and T for J = and J = 0.4|t|, 
respectively (recall that we actually plot (—1) x S*{T) 
to facilitate comparison with NCO). The prediction of 
thermopower enhancement is quite clearly visible. For 
dopings X > 0.5, instead of monotonically decreasing 
from its upper MH limit the thermopower grows to a 
maximum at approximately l\t\ < T < 2\t\ before be- 
ing pinned by its T = constraint. The enhancement 
is also larger for larger dopings. Again the value of J 
has little effect other than at extremely low temperatures 
which most assuredly suffer from finite size effects. As 
expected the thermopower obtains its MH limit when T 
reaches 5\t\ < T < 6\t\ just as for the positive hopping 
situation. 



The origin of the thermopower enhancement stems 
from the transport term Sfj.{T) as expected. Fig. |6^- 
b shows the Mott-Heikes formula for the thermopower 
along with the full S*{T) for the negative hopping case. 
The MH term clearly underestimates the magnitude of 
the enhancement and nearly misses it all together. Ob- 
viously the transport term has a much more sensitive 
dependence on the lattice topology than the MH term. 

The magnitude of the enhancement is also quite strik- 
ing. A value of thermopower greater than 150 iiV/K is 
already anomalously large for a seemingly metallic sys- 
tem such as this. However, upon switching the sign 
of the hopping parameter a value of the thermopower 
of nearly 350 ijlV/K is obtained. Clearly, one should 
consider the intermediate temperature transport effects 
when attempting to discover / design la rge thermopower 
materials; for NCO (|i| - lOOi^ ^ ' ^ ' ^^ ' ^ ^ the intermediate 
temperature range corresponds to approximately room 
temperature. 

The intermediate temperature enhancement of the 
thermopower is very similar in magnitude and shape to 
the enhanced thermopower of NCO at high dopings ob- 
served recently by Lee, et. Although there is, at 
present, no reason to believe that the experimental sys- 
tem has an inversion of the sign of the hopping as the 
doping is increased from the Curie- Weiss metallic phase, 
the similarity between our calculation and the data is 
striking. 

For completeness. Fig. [7|a-f) displays the full fre- 
quency dependence of S{uj^ 1) compared to S* (T) for the 
negative hopping. The same general behavior is shown 
compared to the positive hopping situation in Fig. Qa- 
f). The main difference between the two cases is that the 
deviation of S{uj,T) and S*{T) is slightly more extreme. 
However, similarly, this is seen to occur at low tempera- 
tures, frequencies, and dopings, and is most likely a finite 
size artifact. 



V. LORENZ NUMBER 

The Lorcnz number is an important thermoelectric 
quantity, not least of which, is its importance in the di- 
mensionless figure of merit (FOM). The FOM is an ex- 
tremely important quantity when determining the tech- 
nological usefulness and performance of a thermoelec- 
tric material where values of FOM in excess of unity are 
highly desirable. 

Recall the definition of the Lorenz number from Eq. [3] 



Taiio,T) 
which can be simplified as 



a{u;,T)j 



(23) 



(24) 



Ta{u;,T) 

with k{uj,T) = k{uj,T) - {2fi{T)/qe)^{ij,T) + 
{n{T)/VTqefa{uj,T) being used to define k(w,T). The 
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x=0.83. t<0, J=Oltl 



x=0.75, t<0, J^OItl 



(S((0)-S )[HV/K 




(S((0)-S )[HV/K] 

2 




(a) 



(b) 



x=0.67, t<0, J^OItl 



x=0.58, t<0, J^OItl 



(S((0)-S )[HV/K] 




(S((0)-S )[HV/K] 

2 




(c) 



(d) 



x=0.17, t<0, J=Oltl 



x=0.083,t<0, J=Oltl 



(S(m)-S )[|1V/K] 



(S(CO)-S )[HV/K] 





(e) 



(f) 



FIG. 7: (color online) S{lj, T) — S*{T) as a function of frequency ui and temperature T for J — 0, negative sign on the hopping 
t < (corresponding to a fiduciary hole doped C0O2 compound) and dopings (a) x = 0.83, (b) x = 0.75, (c) x = 0.67, (d) 
X = 0.58, (e) X = 0.17, and (f) x = 0.083. The general frequency dependence is quite similar, albeit larger (discussed further 
in the text) to that shown in Fig. |4] for positive hopping. 



chemical potential drops out of this formula entirely, 
hence, one could calculate L(u),T) in a purely canon- 
ical ensemble. For non-interacting electrons it is easy 
to show that at T = the Lorenz number is equal to 
Lq = (nkB/V^qe)^ ■ This result is simply the familiar 
Wiedemann-Franz law, c.f. Ref. I8|29l The way this 
number is obtained in the non-interacting case via our 
formalism is that at T = there is a delicate balance 



similar to that in the thermopower, namely. 



lim 



T 



-IT 



a(u;,0). 



= 



(25) 



The value of comes from the temperature dependence 
of both terms. As T ^ both k(w, T) /Ta(u}, T) and 
(7(ti;, r)/(T(a;, T))^ behave quadratically in temperature 
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quite satisfactorily determine the accuracy of using the 
high frequency expansion of L*{T) in place of L{uj,T). 
Secondly, we can determine generally qualitative aspects 
of the Lorenz number for strongly correlated systems as 
both a function of temperature and doping. Our quali- 
tative determination of L{uj,T) also allows us to look at 
the dimensionless figure of merit which we investigate in 
section |VI| again at a qualitative level. 

The high frequency expansion of L(uj, T) is 



L*{T) 



T2I 



(27) 



(a) 




(b) 



FIG. 8: L*{T)/Lq as a function of doping x and tempera- 
ture T for positive hopping t > (corresponding to NCO af- 
ter a particle-hole transformation) where Lo — {i^kB/VSqe)'^ ■ 
Panel (a) and (b) are for J = and 0.2|t|, respectively. 



and the difference in the coefficients multiplying each 
quadratic term is equal precisely to Lo- 
in our finite sized system we have neither the delicate 
balance of the T = behavior nor the quadratic low 
temperature behavior present. Using the same "trick" 
as we used for the thermopower we force each term to 
separately vanish at T = 0, i.e., 



LiLo,T) 



1 l Tk{uj,T) 

r2 \ a(c.,r) 



Tk{uj,0) 



a{u;,0) 

1 I a(^,r) 



(26) 



to keep L{u),T) from badly diverging. At low temper- 
atures the exponential behavior caused by the discrete 
nature of the energy levels of the finite system rears its 
head causing problems as T ^ 0. 

At present we are only able to consider a couple of is- 
sues regarding the Lorenz number. Firstly, we can still 



where similarly to T) we have defined Qxx through 

Qxx = Qxx - (2/i(T)/g,)$:,^ + {fi{T)/qe)^Txx. After forc- 
ing each term to vanish at T = Eq. |27] becomes 



L*{T) 



1 (e^^(r)) (9^^(o)) 

T2 ] (tUT)) (t,,(0)) ^ 

2^2 



.my- 



{TxxiO))' 



(28) 



In Figs.jS^L-b and[9|L-f we report results for L*{T)/Lq, 
and (L(a;,T) — L*{T))/Lq for the case of the hopping 
relevant to NCO. The dotted black line in Fig. [8^-b is 
to indicate unity. It is a reasonable assumption that at 
T — Q a. thermodynamically large system would have 
a finite value of L*{T) at T = 0. Whether this finite 
value is equal to is an open question and one we are 
unfortunately not able to shed light upon at this time. 

For the non-interacting case the behavior of L*{T) 
quickly deviates from Lq where it is pinned at T = 
as a function of temperature. For the interacting case 
shown here, the Lorenz number quickly decays to very 
small values as T increases as well. However, our results 
indicate that the intermediate temperature behavior of 
L*{T) generally grows with decreasing doping x, i.e., in- 
teractions evidently increase the Lorenz number. The 
effect of J is much harder to discern due to presumably 
finite size effects and low temperature divergences. J 
seems to have little effect except for the highest dopings 
calculated x =0.17, and 0.083. 

Fig. [9] is similar to Fig. |4] and [7] For dopings x > 0.58 
the frequency dependence of L{ijj,T) is very weak and 
only on the order of approximately ^ (2 — 3)% or less. 
At the low dopings of x =0.17, and 0.083 there is again a 
much a stronger frequency dependence that is likely due 
to finite size effects and not intrinsic to the t-J model. 

With confidence in the weak frequency dependence of 
the thermopower and the Lorenz number we proceed to 
calculate the dimensionless FOM for the situation appli- 
cable to NCO. 
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x=0.83, t>0, J=0.2ltl 
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(b) 



x=0.67, t>0, J=0.2ltl 



x=0.58, t>0, J=0.2ltl 
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(c) 



(d) 



x=0.17, t>0, J=0.2ltl 



x=0.083, t>0, J=0.2ltl 




(L((0)-L )/Lo 



>10 




(e) 



(f) 



FIG. 9: (color online) {L{lu,T) — L*(T))/Lo as a function of frequency uj and temperature T for J = 0.2jtj, positive hopping 
t > (corresponding to NCO after a particle-hole transformation), and dopings (a) x — 0.83, (b) x — 0.75, (c) x = 0.67, (d) 
X = 0.58, (e) X — 0.17, and (f) x = 0.083. The frequency dependence is weak for dopings x > 0.58 ((a)-(d)). For x =0.17, and 
0.083 ((e) and (f)) there is strong dependence occurring at small u and T. Note that in (e) and ((f) we are only considering 
temperatures T > 0.5\t\ as the Lorenz number badly diverges below that temperature. This is most likely a finite size effect 
and not a intrinsic property of the t-J model at these dopings. 



VI. FIGURE OF MERIT 

A value of the dimensionless figure of merit greater 
than one is indicative of a good thermoelectric material 
and is therefore desired. Of course, throughout this work 
we have neglected the lattice contribution to the Lorenz 



number which will add to our calculation of L{ui, T) and 
consequently serve to decrease the FOM that we calculate 
below. 

Recall that the FOM is given in Eg. [4| while the high 
frequency expansion is given by Eq. |16[ The numerator 
is just the square of the thermopower and vanishes at 
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the T = 10|<:| plane is the "normalized" square of the 
MH term, i.e., (ln(2a;/(l — x)))^ indicating an interesting 
doping behavior of the the relative magnitudes of the 
figure of merit. This behavior is most certainly due to the 
MH term of the thermopower. Nonzero J has a nearly 
negligible effect accept at the lowest dopings where it 
serves to slightly reduce the highest temperature value 
of the FOM. 

The full frequency dependence of Z{ijj,T)T is under- 
standably weak (not shown) since both the thermopower 
and Lorenz number have been shown to be weakly fre- 
quency dependent. 



(a) 



VII. CONCLUSION 








(b) 



FIG. 10: (color online) Z*{T)T as a function of doping x 
and temperature T for (a) J = and (b) 0.2|f|. The red 
squares indicate the point at which Z*{T)T = 1 while the 
black dotted line indicates unity as well. 



zero temperature and eventually obtains the square of 
the MH limit, i.e., [{%&^iV/K)\n{2x/{l-x))f. However, 
the denominator is the Lorenz number which starts out 
finite at T = and quickly decays to zero as T — )■ oo. 
Therefore we expect the FOM to begin at zero and grow 
without bound as T increases. 

Fig. [lOk-b shows the FOM for J = and J = 0.2|t| for 



the positive case of the hopping i as a function of both 
the doping and temperature. True to expectations the 
FOM quickly grows from zero at T = to well above 
unity for dopings x > 0.58. At intermediate dopings 
(0.5 > a; > 0.17) the FOM remains very small up to the 
largest temperatures calculated and never reaches unity. 
This is understandable when one considers that the ther- 
mopower in this regime has a small absolute value for 
all temperatures and the Lorenz number only serves to 
diminish this value. At the smallest doping calculated 
{x = 0.083) the FOM just reaches unity as T - 10|i|. 

Note that in Fig. [lO] we have indicated the point at 
which Z*{T)T equals one by a red box. Projected onto 



In this work we have established the general validity of 
the high frequency expansion of the thermopower, Lorenz 
number, and figure of merit for strongly correlated elec- 
tron models. This high frequency expansion is much 
simpler to consider than the full Kubo formalism and 
yet complicated enough to capture the full interaction 
effects. This established validity should provide a bench- 
mark encouraging the theoretical community to obtain 
useful approximate methods to calculate these high fre- 
quency formulas. 

We also provide theoretical evidence supporting the 
authors' previous calculations in re gard to the Curie- 
Weiss metallic phase of NCOi^liil being qualitatively 
and quantitatively described by a two-dimensional t-J 
model on a triangular lattice. Further, we have provided 
new predictions for sodium cobalt oxide concerning the 
Lorenz number and figure of mer it (Sees, fvl and [Vl|. 
In Ref. I13|14l the present authors used an experimen- 
tal hopping parameter of |t| ~ lOOK taken from both 

tpm 



photoemission experiments 



and chemical potential 



measurements^ indicating NCO to have a very narrow 
band system. Therefore, one should scale the temper- 
ature in this work by this value of the hopping when 
making experimental comparisons. 

The predicted thermopower enhancement discussed in 
Sec. |IV 2] should stimulate an experimental search of lat- 
tice based strongly correlated materials to provide useful 
thermolectric materials. Perhaps it is possible to cus- 
tom design a high thermopower material armed with the 
knowledge of the behavior of the transport term alone 
for strongly correlated systems. 
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APPENDIX: OPERATOR EXPRESSIONS 

This appendix presents formulas taken from Ref. lll2l3l for completeness and clarity and additional ones used in the 
present work. 

In discussing the explicit forms of the operators we introduce the so-called Hubbard operators to act to simplify 
the formulas somewhat. The Hubbard operators are defined as follows: 

Xaoirl ^ 4 (A.l) 
Xo^if) = apa (A.2) 

X,,,ir} = 45,-,, =X„o(r0^o.'(O (A.3) 

Xoo(rO = a.>4 - Xo^ir)X,o{r^ . (A.4) 

They have modified anti-commutation relations 

{4'CrV'} = {X^oir^^XoAr")} 

= SpP'{X„„,{r') + 5,„,Xoo{7^)) 

= 5ppY,,,{r') (A.5) 



where the last line defines Ya'cr{'r). Using this notation the charge and energy currents defined in Eq. [sjand Eq. 10 
are given by 

Jx = - lim -^[K{kx),qeh{~kj:)] ^ iqet'^r]xX^o{r + rf)Xo^{f) , (A. 6) 



— ^0 dk'T 

rrj(T 



and 



if = - lim A[f(fc,) if(-fc.)-h;7(-M] 

k:^^o dkx 2 

= — ^ X! i'nx + 'n'x)Yaa'{r + 7l)X„o{r + fj+ff)Xoa,{T^ + —^{i^^fl{r + 7^,f^ 

rr)r}' acr' rrjrj'a 

+ [{Vx + 2ri'JSp+,j+ij, + {f]^ - 2r]'JSp+,f] ■ fl{f+ rf, f)} , (A.7) 

respectively. 

The stress tensor is simply 

Txx = - lini -l^[Jx{kx),qen{-ka:)] = qlt'^rilXcyo{r + fi)Xocr{r) , (A.8) 

kj^~^0 (JjlVrj' 

Tr\(7 

while the modified thermoelectric, and thermal operators {^xx t^^d ^xx) have the form 

^xx = -1™ -^[Jx{kx):H{-kx)] 
K.J. *o ditZx 

= ~9e^ X! i'nx + V^fY„„,{r + rii)X„Q{r + fj+Tf)Xo„:{T^ + -^^{i^^[i^^il{^ 

rr}rj' oa' rip]' a 

+ [{Vx + 27y^)5,-+,y+,7 -f (77, - 2r^'^)Sp+.^,] ■ pi{r + rj, r)} , (A.9) 

and 

kxx = - Vim^^[J^{kx).H{-kx)] 
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rr)r}'r]" o 

-{X^o(r + r/)Xos(f + ?7 + rf' ) ~ /i.c.}X„o(r + ^')^0a(r1] 

+ E + + ^^')(^- + ^'^x + 0^<.a'(r + ^r,,v'(r + 77 + rr)X,,.o(r + ^+ rf + ^')^o<x'(r)} 
+ 1^ E fe + »?x){E[('?^-^x)^ao(r + ^+jf)Xo<.'(r)[X<,.s(f+^)Xs^(f+^+?f') 

Tr\r\' r\" era' 

-Xs'a{r + rf)X^,s,{r + r^+rf')] - (77^ + 77^)y^^.(f + ?7')X3o(r + + ^)Xoa'(r)^a-?(r + + + ?/') 
-(ry, + ii',)Y,a'{r + rf)X,n{r + + ^)Xo^- rf') 

H^x -ri'^ + 2ril)[X^^{r + 77 + if')X^,-^{r + ff) - X^,^{f+ f] + ff')X^>^{f+ ff)]X^o{r + ff+ rf)Xo„,{r) 

-{Vx + + '^v'^)X^a(r + ri+ rf + rf')Y„„,(r + 'Fi)X^o{f+ q + rf)Xa„,{r) 
-iVx + v'x'- ^Vx)Xa'a'{r + rf')Y,^'{r + rf)X,o{T^+ rf+ 7f)Xas'{r) 

-{Vx +'nx'- Vx)Yaa'{r + ff)X^o{f+ f] + ff) Xq^' {r)[np+fj" a' " np+ff"a'] 

-{Vx + Vx + Vx)^oa'{r + ff)[np+ff+fj>+fj"a - np+,j+ff>+f}"g]Xaa{f+ rf + rf)Xo„>{r)] 
+{2t]^ - 2t]'^ + 2t]'^) ^[nf+ff+fj„a - np+fi+fj>,s]Xg„{f+ ■q)X„o{r + rj + Tf)Xng{)} 



tJ 

E {E^^«/37[^xM"(^+^,r=){5^+^„ 

rrjrj'rf" afS'y 



(77, - 277^ 

4 '^r+rf" 4 "^f+^+rj" 2 '^f+fi'+fi" 



Vx Vx QJ 



^ r+rj" ^ f+rj+T?" 



f+ff+ff'+if'' 



} 



+Vx 
+Vx 



iVx - 20 + 277^0 00 (7?. - 27?^ - 277^0 CO 

^ ^f+rj" ^ ^r+ff+rj" 2 "^r+jj'+rf" 

(^x - 277^0 (^x + 277^0 oa (77^ + 277^ + 277^0 



/x^(r + 7?,f)52;^-, 



M'^(f+,?,r-)52+^^+^. 

+(772: + 2r]'^)[Spj^fi_^_fj, {f + f],f) Sj_^_fi_^_fj„ - ^Sjj^f^,, - — 



+ 



(77, + 277^') _ (r?, - 27?;') ^„ (77, + 277; + 277; 

' '-'r+ff+if" 



+ 



^ - r-i-Tj-l-T/ ■ -r+T?" 2 

% c7 _ »fe (.7 (77x - 277; 

(77x + 27?;0 (??x - 20 , - 27?; - 27?^') ca 

^ '-'r+ff+77" ^ '-'r+ff 



+(77. - 277;')[5.V/(r + 77, r) <| ^SJ^,^.,, - '^S}^,., + 



x) 07 \ 



xJ na Vix "'ixJ qa 1 \'ix —ix —ixJ ga I c/3 in(r-i^f;r\ 



VxiVx-2r]'^) T]x{Vx + 2 77^0 



'^fM-i7+^"| ^0-0(7^+ v)^Oa{r) {52^^, + -S'2^^_,_^,} 



+[(77. - 277;')52^^^" + (r?x + 2rj':)S^^^+f,„]S^+ff+f,,X,o{f+ ff}Xo.{f)} 
iVx - 277^) 



+ 



-{Sp_^_fj,X^o{r + v)Xo^{T^[rj^Sp+fj„ + -nxSp+fj+fj,, 



+ [fe - 27j':)Sp^^„ + (r?, + 2r?;')52^^^+,-"]52+^?'^.o(r + 77')^Oa(r1}]} , 
respectively. In the above expression for ^^x smd Q^x we have used the bond spin operators defined as 



(A.10) 



H^ir,/) = -[X^o{r)Xo^{r')-X^oir^Xoii/)] 



(A.11) 
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Further note that for Eq. A. 10 we have used the definition 

d n 



J^{k,) = - lim -p- \ ^[T(g, + fc,),T(-g,)] + [T{q, + k,),V{^q,)] 
allowing the establishment of the identity (using the Jacobi identity) 



(A.12) 
(A. 13) 

(A. 14) 
(A. 15) 



simplifying the calculation of Eq. |A.10| tremendously. 
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